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Abstract 

The lightest neutral scalar in the Inert Higgs Doublet Model is a natural candidate for WIMP 
dark matter. In this paper, we analyzed the dark matter relic density in the Inert Higgs Doublet 
model. Various theoretical and experimental constraints are taken into account. We found that 
there are five distinctive regions that could provide the right amount of the relic density in the 
Universe. Four out of those five regions have a light particle spectrum which could be studied at 
the Large Hadron Collider. 
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I. INTRODUCTION 



About 20% of the Universe is made of cold dark matter. The origin of the dark matter, 
however, still remains a mystery. While the Standard Model (SM) of particle physics 
has been very successful in explaining the data from almost all of the particle physics 
experiments to date, none of the SM particles can be a good candidate for the dark matter. 
Its existence provides unambiguous evidence for new physics beyond the Standard Model. 

A Weakly Interacting Massive Particle (WIMP) is a promising candidate for the dark 
matter, given that the WIMP relic density is naturally around the observed value [l[ 

fi/i^ = 0.112 ±0.009, (1) 

for a WIMP mass around the electroweak scale ~ 100 GeV. WIMPs also appear naturally 
in many beyond the Standard Model scenarios. One popular candidate is the lightest 
neutralino in supersymmetric models, which provides an example of spin- 1/2 WIMP dark 
matter. Spin-1 WIMP dark matter has been studied in the framework of extra dimension 
models, for example, the lightest Kaluza-Klein photon in Universal Extra Dimension 
models [2|. 

Spin-0 dark matter has been studied 0, 0, [El in the framework of the SM plus an extra 
scalar. The minimal model includes an extra real scalar gauge singlet S that is charged 
under a Z2 symmetry. The only renormalizable couplings are the quartic self-coupling 
XsS"^ and the Higgs-S* coupling XS^WH. It was shown that for a SM Higgs with mass in 
the range of 100 GeV < irih < 200 GeV and dark matter with mass in the range of 10 GeV 
< ms < 100 GeV, viable regions of parameter space exists for the dark matter candidate 
S to provide the right amount of relic density in the Universe jsf. 

Unlike the WIMP, the mass and the couplings of such a scalar dark matter particle 
often needs to be fine tuned to be consistent with the observed dark matter relic density. 
WIMP-type scalar dark matter could appear in a simple extension of the SM with the 
addition of a Higgs SU(2)i doublet. Unlike the SM Higgs doublet, which couples to both 
fermion and gauge sectors, such an extra scalar doublet is inert in the sense that it couples 
to gauge bosons only. Such inertness can be guaranteed by imposing a discrete Z2 parity. 
The lightest inert particle (LIP) is therefore stable. If it is neutral, it can be a good WIMP 
dark matter candidate. 

The Inert Higgs Doublet Model (IHDM)^ was first proposed in the late 70 's [6]. It has 
received recent attention [3| since it could be used to solve the naturalness problem in the 
Standard Model. With a splitting in the masses of the neutral and charged components of 
the inert Higgs fields, the resulting positive contribution to the oblique T parameter allows 
this model to accommodate a SM Higgs with a much heavier mass. Such an inert Higgs 
sector also appears in the recently proposed left-right Twin Higgs model In addition, 
an inert Higgs doublet (or doublets) could be used to explain the small Majorana neutrino 



Strictly speaking, the extra SU(2)l scalar doublet is not a Higgs doublet since it does not obtain a 
vacuum expectation value. Such model is also referred to as "Inert Scalar Doublet Model" in the 
literature. 
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mass via the one-loop radiative seesaw mechanism [ol, electroweak symmetry breaking T3 
grand unification 11 1 and leptogenesis [l^. 

There have been some studies on the dark matter candidate in the IHDM. Ref. 
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studied the dark matter relic density in the IHDM in certain regions of parameter space. 
The right amount of relic density could be obtained for dark matter with mass around 40 
GeV — 80 GeV or larger than 600 GeV. Ref. IJ, ll5| studied the neutrino signatures from 
dark matter annihilation. The continuous gamma ray spectrum from fragmentation and 
monochromatic gamma ray lines were studied in Ref. 13(] and IGj respectively. Positron 
and antiproton signatures were studied in Ref. [l7|. Direct detection has been studied in 



0, [H, 18]. There is also a collider analysis based on the LEP II limit [l^ as well as collider 
signatures of 5*^4 associated production with A — > Sl~^l~ at the Large Hadron Collider 
(LHC) t20|]. 

In this work, we performed a complete analysis of the dark matter relic density 
in the IHDM over the whole parameter space, taking into account various theoretical 
and experimental constraints. The latest results of the collider constraints based on 
sup ersy mmetric process e~^e~ —>■ X1X2 search at the LEP II are imposed [l9|. Unlike 
Ref. j_13| , in which only a low SM Higgs mass rrih = 120 GeV and 200 GeV are considered, 
our analysis also includes studies with a high SM Higgs mass rrih = 500 GeV. In Ref. , 
the mass splittings between if^, A and the dark matter candidate S, {mH± — "^5, "^a— "^5)5 
are fixed to be (50, 10) GeV and (10, 5) GeV for the low and high dark matter mass regions, 
respectively. We study cases when the mass sphttings between if^, A and S are small, in 
which coannihilations play an important role, as well as cases when the mass splittings are 



large. In regions that overlap with those analyzed in Ref. [13(], our results agree with the 



literature. We identify additional regions of parameter space, in which the dark matter 
relic density is also consistent with the WMAP result but were overlooked previously. 
We also present our results in the parameter spaces of physical Higgs masses and Higgs 
couplings, which are convenient to use for studies of collider phenomenology and dark 
matter detections. 

The rest of the paper is organized as follows. In Sec. [Tll we briefiy present the IHDM. In 
Sec. Illlt we discuss the theoretical and experimental constraints on the model parameter 
space. In Sec. HVl we present our results of the relic density analysis. In Sec. |Vl we 
conclude. 



II. THE INERT HIGGS DOUBLET MODEL 

The IHDM is an extension of the Higgs sector of the SM. Besides the usual Higgs doublet 
Hi = HsM, an additional Higgs doublet H2 is introduced: 

( (5+^^)772 )' 

which is charged under SU(2)i x U(l)y as (2,1/2). Unlike the SM Higgs boson, which 
couples to both SM gauge bosons and matter fermions, the extra Higgs doublet H2 couples 
to the gauge sector only. Such couphngs can be guaranteed by imposing a Z2 symmetry 
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(sometimes also called matter parity) under which all particles except H2 are even. In 
particular, under the Z2 symmetry: 

Hi — > Hi, H2 —H2. (3) 

While Hi obtains a vacuum expectation value (VEV) v/y/2 = 174 GeV as in the SM, 
H2 does not obtain a VEV: {H2) = 0. The Z2 symmetry is, therefore, not spontaneously 
broken. The lightest particle in H2 is stable and could be a good dark matter candidate. 

The most general CP-conserving potential in the Higgs sector that respects the Z2 
symmetry can be written as 



V = i4\Hi\^+i^l\H2\^ + Xi\Hi\^ + X2\H2\^ + X3\Hi\^\H2\'' + Xi\HlH2f + 



^{HlH2r + h.c. 



(4) 

Notice that the usual mixing term /if2-ffi-ff2 is forbidden by the Z2 symmetry. After 
electroweak symmetry breaking, three degrees of freedoms in Hi are eaten by massive 
gauge bosons and Z. We are left with one physical Higgs boson h, which resembles 
the SM Higgs boson, as well as four inert scalars: the CP even one S, the CP odd one A 
and a pair of charged ones H^. The mass of h is related to Ai via 

ml = -2fxl = 2XiV^. (5) 

The masses of S, A and H^ are related to parameters in the Higgs potential as 

ml± = + Xsvy2, (6) 
ml = fil + (A3 + A4 + A5)t;V2, (7) 
ml = + (A3 + A4 - X5y/2. (8) 

(9) 

We define the mass differences 61 and 62 as 

iX4 + X5)v^ X^v^ 
di = mH±-ms = -— ■ r, d2=mA-ms = --, ■ r. (10) 

2[mH± + ms) [mA + ms) 

It is obvious that A4 and A5 control the mass splitting between the charged and neutral 
CP even states, while A5 also controls the mass sphtting between the CP odd and CP even 
states. In our analysis below, we assume that the CP even scalar S is the LIP, therefore the 
dark matter candidate, giving 5i 2 > 0. The numerical results of the relic density analysis 
are similar if A is the LIP dark matter. 

The Higgs potential in Eq. (j4]) has seven parameters: 

(/^15 /^2' -^1; -^2, A3, A4, A5). (11) 

They could be replaced by the Higgs VEV f , physical Higgs masses, mass splittings, and 
a sum of quartic couplings Al = A3 + A4 + A5 as 

(i;,mft,m5,5i,52, A2, Ai). (12) 
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In particular, shows up in the couphngs of SSh and SShh, which is relevant for dark 
matter annihilation. It is therefore convenient to pick as a model parameter. The 
quartic coupling A2 only shows up in self-couplings involving 5", A and H"^. It does not 
play an important role for the dark matter analysis that we present below. 



III. THEORETICAL AND EXPERIMENTAL CONSTRAINTS 

There are various experimental constraints on the IHDM from direct collider searches, 
indirect electroweak precision test and dark matter direct detections. 

• W and Z decay widths 

Light H^, S and A could lead to deviations of Z and W decay widths from the SM 
value. On the other hand, T^z have been measured precisely at the LEP and the 



Tevatron [21[, which agree well with SM predictions. Therefore regions in which the 
decay processes —>■ SH^/AH^ and Z H^H^ / SA are kinematically allowed 
have already been excluded. These translate into constraints for ms and 61^2 as 

2ms + 5i > mw, 2ms + Si + 62 > mw, 

2ms + 2(5i > mz, 2ms + ^2 > mz- (13) 
Direct collider searches 

Light neutral and charged Higgses have been searched for at the LEP and the 
Tevatron. Limits from conventional searches for if^, S and A, however, do not 
apply for Higgses in the IHDM since those searches rely on decays of Higgses into 
fermion pairs and/or Higgs production via top quark decay. Neutral and charged 
Higgses in the IHDM, on the other hand, do not couple to fermions. In particular, 



the charged Higgs has been searched for at the LEP and the Tevatron [22, l23[ . A 



lower mass bound of 74 — 79 GeV at 95% C.L. is obtained at the LEP considering 



cs, r+z/. A more recent search at CDF [23| studied the charged Higgs produced 
in the top quark decay t H^b, with further decaying into a pair of quarks, 
leptons, or W~^(f). The bounds on the charged Higgs mass from those searches do not 
apply to Higgses in the IHDM since Higgs- fermion-fermion couplings are absent. 

Experimental signatures for S, A and H^, however, is similar to those of neutralinos 
and charginos in the Minimal Supersymmetric Standard Model (MSSM), as S 
appears as missing energy at colliders, similar to the lightest supersymmetric 
particle(LSP) in the MSSM. In particular, searches of e+e- ^ x?X2 at the LEP 
II can be interpreted as searches for e~^e~ — > 5*^4. Therefore, the null result for 
neutralino and chargino searches at the LEP II can be used to set limits on m^^, 
mA and ms in the IHDM. A recent analysis [l^ argued that a direct application of 
the upper limit on the LEP II X1X2 cross section to the IHDM is oversimplified due 
to the difference between the MSSM process X1X2 and IHDM process 5*^4. Based on 
DELPHI analyses of e+e~ —>■ X1X2 with X2 ~^ XiQQi A*'*'/^") e"'"e~, Ref. [ll] determined 
the efficiencies for the corresponding MSSM and IHDM processes after cuts using 
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Monte-Carlo simulations. The ratio of the efficiencies is then used to rescale the 
MSSM cross section upper limit and applied to the IHDM. It is shown that regions 
satisfying ms ^ 80 GeV, ^100 GeV and ^2 > 8 GeV are excluded by LEP 
II MSSM searches. This limit is stronger than previous estimations by the direct 
application of the X1X2 cross section upper limit to the IHDM 0, Hi- ^2 < 8 
GeV, however, limits on scalar masses are much weaker, due to the small mass 
splitting and the resulting soft jets and leptons in the final states. For ^2 < 8 GeV, 
only the LEP I limit, ms + ttia > mz, applies. 



Similarly, searches of the supersymmetric channel e+e^ XiXi the LEP II [25 
can be used to set a bound on mH± . Taking into account the cross section difference 
of scalars versus fermions, a limit of mH± ^ 70 — 90 GeV can be derived from the 
LEP II chargino searches [23 |. 

Electroweak precision test (EWPT) 

Electroweak precision measurements provide strong constraints for any new physics 
beyond the SM. In the global electroweak precision fit to the SM, a light Higgs 
is preferred: rrih = 90tlj GeV, with ruh < 163 GeV at 95% C.L. Oblique 
parameters [23] S, T and U are often used to parameterize radiative correction to 
gauge boson propagators. A heavier SM Higgs contributes positively to 5* while 
negatively to T. Extra Higgs bosons in the IHDM, namely H^, S and A, also 
contribute to the S and T parameters, which could possibly cancel the effect of a 
heavy SM Higgs boson [3]. Therefore, a heavy SM Higgs can be accommodated 
in the IHDM. In our analysis, we require the overall contribution to the S and T 
parameters from extra scalars in the IHDM to fall within the 68% C.L. ellipse in S*— T 
plane. Since we considered possible large mass splittings 61^2, where the approximate 
formulas for AS* and AT in Ref. Q might not be valid, we used the full expressions 
for the contributions to AS" and AT for a Higgs doublet 0| when imposing this 
constraint. 

Fig. [T] shows the allowed region for mass splittings 5i and 82 given the constraints 
on the oblique parameters 5* and T from electroweak precision measurements, for 
several different choices of ms- The dependence on is small, especially for large 
6i^2- For mh = 120 GeV (left plot), the allowed region falls along the diagonal 
direction of 61 ~ 62- The mass splittings 61^2 could both be large as long as these 
two mass splittings are close to each other. For mh = 500 GeV (right plot), a large 
contribution from the inert Higgs doublet is needed to cancel the large positive S 
and negative T contribution from a heavy SM Higgs boson. Therefore, the value for 
61 is constrained to be quite large 61 ^ 150 GeV, while the value for 62 could still be 
as small as 0. 

Dark matter detection 

Dark matter direct detection excluded spin-independent dark matter— nuclei scatter- 



ing cross section up to about 10~'^^cm^ at 90% C.L. |28|, |29|, which is 7—8 orders of 
magnitude larger than dark matter — nucleon scattering via Z exchange. This limit, 
however, can be easily avoided by introducing a small mass splitting between S and 
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FIG. 1: Plot of the 68% C.L. allowed region in mass splittings 6i and 82 given the constraints on 
S and T from precision electroweak measurements. The value for the SM Higgs mass is set to be 
120 GeV for the left plot and 500 GeV for the right plot, ms is taken to be 10 GeV (solid curve), 
40 GeV (dashed curve), 75 GeV (dash-dotted curve) and 600 GeV (dotted curve). 



A: ^2 ~ a few hundred KeV, given the typical kinetic energy of the dark matter and 
the momentum transfer between the dark matter and the scattering nuclei in such 
scattering processes. 

The spin-independent dark matter— nucleon scattering cross section via SM h- 
exchange is given by |7(] 

471 V rns J \ml J 

where m^v is the nucleon mass, is the reduced mass of the dark matter— nucleon 
pair, and / is used to parameterize the nucleonic matrix element. The typical range 
for / is taken to be 0.14 — 0.66 [si. In our analysis, we take / to be 0.14 to be 



on the safe side. Latest results from XENONIO, CDMS, CRESST, CoGeNT and 



TEXONO are used in our analysis [28|, l29|, l31|, |32|, l33|, l3^. Portions of parameter 
space in small ms, rrih and large |A£,| region are excluded by the upper limits on the 
spin-independent dark matter— nucleon cross section. 

The bounds from indirect dark matter detection (gamma rays, for example) is weak 
in the IHDM. Moreover, those bounds are subject to large astrophysical uncertainties 
involved in those observations. Therefore, we don't consider constraints from indirect 
dark matter detections. 

In addition, we impose the following theoretical constraints. 
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• Vacuum stability 

We require the vacuum stability of the Higgs potential at tree level, which leads to 

Ai,2 > 0, (15) 
A3,A3 + A4- IA5I > -2v/A^. (16) 

• Perturbativity 

We require that corrections to the beta function of Ai from non-SM quartic couplings 
A3^4^5 is less than the 50% of the SM term 24Af [7]. This amounts to the constraint: 

A^ + (A3 + A4)' + A^ < 12X1 (17) 

The evolution of the remaining quartic couplings does not lead to extra constraints. 
In addition, we require the quartic coupling A2 to be in the perturbativity region: 

A2 < 1. (18) 



IV. RELIC DENSITY ANALYSIS 



We analyzed the relic density in the IHDM using the program MicrOMEGAs [35[. This 
program solves the Boltzmann equation numerically, using the program CalcHEP [s^ to 
calculate all of the relevant cross sections. When the mass splittings between the dark 
matter candidate and other particles are small, coannihilation effects are also included. 

In our analysis, we fixed the SM Higgs mass to be 120 GeV (preferred by the EWPT in 
the SM) or 500 GeV (preferred by naturalness, and could be consistent with the EWPT 
with large splittings between H^, S and A.). We fixed A2 = 0.1 since it does not enter 
into the dark matter relic density analysis. We studied both the cases of small mass 
splittings and large mass splittings. When 5i{52) is small, coannihilation between S and 



H'^{A) is important. Part of our study overlaps with the analysis in Ref. [13|, which, only 
considered small Higgs masses rrih = 120 GeV and 200 GeV with fixed mass splittings 
(^1, ^2) = (50, 10) GeV for the low mass region and {61, 62) = (10, 5) GeV for the high mass 



region. Results in Ref. [13[ were presented in jj2 — rns plane. We present our results in 
\l — iTLs plane instead, which is more transparent since ms and A^ are the two relevant 
parameters for the dark matter analysis, as well as studies of dark matter direct and indirect 
detections. In particular, the value of A^ is closely related to the neutrino and gamma 
ray flux from dark matter annihilation, and the cross section for dark matter— nucleon 
scattering. Larger typically leads to enhanced flux and increased dark matter— nucleon 
scattering cross sections, and the inert Higgs dark matter has a better chance to be detected 
at future dark matter detection experiments. 

We discuss below in detail two mass regions of that could provide the amount of 
dark matter relic density consistent with the WMAP result at the 3cr level: (A) low mass 
region where ms < 100 GeV, and (B) high mass region where 400 GeV < ms < around a 
TeV. 
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A. Low mass region 
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FIG. 2: Dependence of the relic density on ms for m/i=120 GeV and 6i = 50 GeV. In the left 
plot, 82 is chosen to be 10 GeV (solid line), 8 GeV (dashed line) and 5 GeV (dash-dotted line) 
while Xl is fixed to be 0.01. In the right plot, Xl is chosen to be 0.01 (solid line), 0.1(dashed 
line), 0.2 (dash-dotted line), 0.3 line), 0.5 ("o" line) while 82 is fixed to be 10 GeV. The 
horizontal band indicates the 3a region that is consistent with relic density measurement from 
WMAP: 0.085 < nh'^ < 0.139. 

To illustrate the dependence of the relic density flh"^ on nis, Xl and mass splittings 61^2, 
we present in Fig. [2] curves of Qh^ vs. ms for various choices of 62 (left plot) and (right 
plot). The horizontal band indicates the 3a region that is consistent with relic density 
measurement from WMAP: 0.085 < ilh"^ < 0.139. A light SM Higgs mass nih = 120 GeV 
is used in both plots. 

The red curve in the left plot of Fig. [2] corresponds to (^i, ^2) = (50, 10) GeV. Therefore, 
coannihilation between 5* and A is important while coannihilation between 5* and is 
not. For small ms, SS bb via SM h exchange dominates, with cross section proportional 
to A|. The cross section is typically small due to the small bottom Yukawa coupling, which 
leads to relic density too big that overdoses the Universe. When ms gets larger, SA — > qq 
via Z exchange becomes more and more important. The relic density enters the WMAP 3a 
region for ms around 35 GeV. The coannihilation cross section maximizes at the Z-pole: 
ms + mA ~ mzi corresponding to the dip around m5 ~ 40 GeV. The coannihilation cross 
section decreases when ms increases away from the Z-pole region, which makes the relic 
density falls back to the allowed region. As ms gets larger, SS bb annihilation via h 
starts to dominate and the relic density enters the 3a region again. For ms ~ mh/2, SS 
annihilation hits the /i-pole, corresponding to the second dip around ms ~ 60 GeV. The 
annihilation cross section gets smaller once ms leaves the /i-pole region. When ms ^ 70 
GeV, SS — >■ WW dominates. The annihilation cross section quickly increases and the relic 
density drops below the WMAP observed value. 
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The dashed and dot-dashed curves in the left plot of Fig. [2] shows the relic density 
dependence for 62=8 GeV and 5 GeV respectively. Coannihilation effects get stronger for 
smaller mass splittings. Therefore, for most of the ms region between 40 — 60 GeV, the 
coannihilation cross section is too large and the relic density is too small. 

Curves in the right plot of Fig. [2] correspond to Al=0.01, 0.1, 0.2, 0.3 and 0.5, 
respectively, while (61,62) is fixed to be (50,10) GeV. Similar Z-pole and h-pole features 
appear. The relic density is smaller for larger A^, since SS annihilation via /i-exchange is 
increased due to the increased SSh couphng. 




FIG. 3: WMAP 3a allowed region (enclosed by blue curves) in ms — \l plane for mh=\2Q GeV. 
The mass splittings are chosen to be {61,82) = (50,10) GeV (left plot) and (50,8) GeV (right 
plot). Shaded regions are excluded either by LEP I+II searches (yellow, light shade), electroweak 
precision constraints (orange, medium shade), dark matter direct detection (purple, medium-dark 
shade), vacuum stability (red, dark shade along bottom), and perturbativity (hatched region). 

Fig. [3] shows the WMAP 3cr relic density allowed region (enclosed by two blue curves) 
in the nxg - \l plane for nis < 100 GeV with rrih = 120 GeV for {61,62) =(50, 10) GeV 
(left plot) and (50, 8) GeV (right plot). Shaded regions are excluded by various theoretical 
and experimental constraints, as described in Sec. IIIII 

For {61, 62) =(50, 10) GeV (left plot), the gap around ms ~ 40 GeV corresponds to the 
Z-pole. The gap around 7715 ~ 60 GeV corresponds to the h-pole. The LEP constraint 
(yellow, light shade region) is very strong due to the strong constraints on and ms 
when 62 > 8 GeV. The precision electroweak constraints (orange, medium shade region) 
is weak since 61 > 62 is slightly preferred by the fit to the S — T contour. Given all the 
constraints, only a small region around ms ~ 80 GeV survives. The value for for the 
allowed region, however, could be as large as —0.2. Such a large value of Al would be 
important for generating a large signal in the indirect detection of dark matter. 

The LEP constraints on ms and m^, however, are weakened for small mass splitting 
62 ^ 8 GeV. For such a small mass splitting, as low as around 40 GeV is still allowed. 
In the right plot of Fig. [31 the allowed parameter space is given for {61, 62) =(50, 8) GeV. In 
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most of the mg region between 40 GeV and 60 GeV, the rehc density for the dark matter 
is too small due to the large coannihilation SA cross section. However, there is a viable 
region for 60 GeV < < 80 GeV with A^, in the region of ±0.1. 
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FIG. 4: WMAP 3cr allowed region (enclosed by blue curves) in ms — Xl plane for m^=120 GeV. 
The mass splittings are chosen to be {5i,52) = (10,50) GeV (left plot), (10,10) GeV (right plot). 
Shaded regions are excluded by various theoretical and experimental constraints; see caption of 
Fig. [3] for explanation. 



The left plot of Fig. H] shows the allowed relic density region for {5i,52) =(10 , 50) 
GeV and uih = 120 GeV. The dark matter annihilation behaves similarly to the (50,10) 
case, only that the previous SA coannihilation via Z is replaced by SH"^ coannihilation 
via W^. The LEP constraints (yellow, light shade region) becomes weaker comparing 
to {61,62) = (50,10) GeV case, due to the weaker constraints on mj^±. The precision 
electroweak constraints (orange, medium shade region), however, becomes stronger. Given 
all the theoretical and experimental constraints, a region with 60 GeV < < 80 GeV 
and —0.1 < Ai < 0.1 survives. 

When 61 and 62 are both small, the coannihilations between S, A and if ^ are relevant 
and the allowed 3a regions shrink due to the enhanced coannihilation cross sections. For 
(5i,52) = (10, 10) GeV, there is no allowed region that survives given all the theoretical 
and experimental constraints, which is shown in the right plot of Fig. |H For ^2 ~ 8 
GeV, although the LEP II constraints is weaker, the WMAP 3a region for ms > 40 GeV 
completely disappears due to strong coannihilation effects. 

Fig. [5] shows the allowed relic density region for ruh = 120 GeV when both 5i 2 are 
large: {61,62) = (50,50) GeV (left plot) and {61,62) = (70,70) GeV (right plot).' The 
gaps in the low mass region corresponding to the Z-pole or l^-pole disappear, since the 
mass splitting is too large for the coannihilation process to be important. The /i-pole, 
however, still survives due to the low value of rrih = 120 GeV that we pick. The LEP 
I-I-II constraints (yellow, light shade region) are weaker because of the large value of 62- 
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FIG. 5: WMAP 3a allowed region (enclosed by blue curves) in mg — Xi plane for 171^=120 GeV. 
The mass splittings are chosen to be {61,62) = (50,50) GeV (left plot), (70,70) GeV (right plot). 
Shaded regions are excluded by various theoretical and experimental constraints; see caption of 
Fig. [3] for explanation. 



The perturbativity constraints (hatched region), however, are stronger, since larger mass 
splittings 61^2 correspond to larger values for A4,5. 

For splitting (^1,52) = (50,50) GeV, a region with 55 GeV < ms < 90 GeV and 
—0.2 < Ai < is consistent with the WMAP 3a region. When the mass splittings 
get larger, the perturbativity constraint gets stronger (the hatched region shifts to the 
left) while the LEP II constraint gets weaker (the yellow, light shade region also shifts 
to the left). For {61,62) = (70,70) GeV, a region with 30 GeV < < 50 GeV and 
—0.15 < Xl < —0.05 opens up. For {61, 62) = (90, 90) GeV the perturbativity constraint is 
so strong that no mass window survives all of the constraints. 

Fig. [6] shows the allowed relic density region for a high value of the SM Higgs mass 
m/j = 500 GeV. Electroweak precision constraints require a large value for 61 (as shown in 
Fig. [1]) in order to compensate for the AS > and AT < contributions from a heavy SM 
Higgs. The left plot of Fig. [6] shows the allowed relic density region for {61, 62) = (250, 110) 
GeV. No pole regions appear due to the large mass splittings and large nih. For nis ^ 75 
GeV, SS —>■ bb dominates, which corresponds to the nearly horizontal band of the 3a 
WMAP region. Once ms ^ 75 GeV, SS — > WW opens up, which leads to the nearly 
vertical band of the 3a region. Dark matter with mass around 75 GeV is allowed given all 
of the constraints. Note that the perturbativity bounds are much weaker due to the large 
SM Higgs mass (therefore, large Ai). On the other hand, to obtain the right relic density, 
large A^ is needed to compensate the suppression of the annihilation cross section by the 
large Higgs mass mt. Although A^ could be as large as 3, the indirect detection in this 
region is not promising due to large suppression of the annihilation cross section by m\. 

The numerical results do not change much for smaller values of 62- A Z-pole 
coannihilation region appears for 62 around 10 GeV. The right plot of Fig. shows the 
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FIG. 6: WMAP 3a allowed region (enclosed by blue curves) in mg — Xi plane for m;j=500 GeV. 
The mass splittings are chosen to be {61,62) = (250, 110) GeV (left plot), (180,8) GeV (right plot). 
Shaded regions are excluded by various theoretical and experimental constraints; see caption of 
Fig. [3] for explanation. 



allowed relic density region for (61,62) = (180,8) GeV. The gap around ms ~ 40 GeV is 
due to the Z-po\e. The LEP II constraints on ms and rriA do not apply, due to the small 
mass splitting. This allows a relatively large ms region around 50 — 80 GeV to remain 
open, with — 1 < Al < 3. This region, however, shrinks for smaller 62 due to the stronger 
coannihilation effects. Therefore, no WMAP 3a region survives for 62 less than about 6 
GeV. 

B. High mass region 

Fig. [7] shows the relic density in the ms — Xl plane for the high mass region: ms > 
400 GeV, with m/i=120 GeV. The LEP search bound is irrelevant now since the ms is 
much larger than the direct search limit. The precision electroweak constraint is weak for 
mfi around 120 GeV as long as 61 and 62 do not differ too much. The dark matter direct 
detection constraint is also weak due to the heavy dark matter mass. For small mass 
splittings , {61, 62) = (1, 1) GeV (left plot in Fig. ED , a region of ms ~ 500 - 700 GeV and 
—0.2 < Al < 0.2 is allowed. The dominating annihilation channels are WW, WZ and ZZ. 
The dependence on A/, is introduced by the SSh coupling. For A^ close to zero, a smaller 
value of ms is preferred to increase the annihilation cross section. When \Xl\ gets larger, 
a larger value of ms is needed for the relic density to fall into the 3a band. The region to 
the right of the 3a band overdose the Universe while the region to the left of the 3a band 
corresponds to the under-abundance region. Regions with large values of \Xl\ are excluded 
due to the perturbativity constraints. 

The annihilation cross section grows for large mass splittings. Therefore, the WMAP 
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FIG. 7: WMAP 3a allowed region (enclosed by blue curves) in mg — Xi plane for 171^=120 GeV. 
The mass splittings are chosen to be {6i, 62) = (1, 1) GeV (left plot), (1, 10) GeV (right plot). Red 
region are excluded by vacuum stability while the hatched region are excluded by perturbativity 
constraints. 



allowed region shifts to larger ms for larger 5i 2, as shown in the right plot of Fig. [7] for 
(61,62) = (1,10) GeV. ms > 650 GeV falls into the WMAP 3a region. The result for 
(^1,^2) = (10, 1) GeV is very similar. The region forbidden by perturbativity constraints 
(hatched region), however, shifts to the left for larger 61^2- Therefore, no allowed region is 
left if at least one of 61^2 ~ 12 GeV. 

For a large SM Higgs mass mh = 500 GeV, a large mass splitting 61 > 150 GeV is 
needed to satisfy the precision electroweak constraints. There is no region in ms — that 
survives after all of the experimental and theoretical constraints are taken into account. 

V. CONCLUSION 

We studied the simple extension of the SM Higgs sector when an extra inert Higgs 
doublet is introduced that couples to the gauge sector only. The lighter of the neutral 
components could be a good WIMP dark matter candidate. We explored the parameter 
spaces of the IHDM, taking into account the relic density constraints from WMAP and 
various theoretical and experimental constraints. Table [T] summarizes five distinctive 
regions that could provide the right amount of cold dark matter in the Universe which 
satisfy all of the constraints. 

In regions (I) — (IV) , the dark matter candidate S, along with at least one of the other 
scalars A, is light. Those particles could be pair produced at the Large Hadron Collider 
(LHC) as SA, AH^, SH^ and H^H^ with cross sections around the fb level. Heavier 
scalars A and could decay into the lightest one S via on-shell(or off-shell) Z*^*) and 
W'^*\ which further decay into quarks, leptons and neutrinos. There are typically large 
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rrih 


ms (GeV) 


Al 


5i, 62 


(I) 


light dark matter 


low nih 


30 - 60 


-0.15 to 


50 GeV< (5i ~ ^2 < 90 GeV 


(11) 






60 - 80 


-0.2 to 0.2 


at least one of ^i, 62 is large 


(III) 




high nifi 


50 - 75 


-1 to 3 


large 61 and small 62 < 8 GeV 


(IV) 






~ 75 


-1 to 3 


large 61 and 62 


(V) 


heavy dark matter 


low nih 


500-1000 


-0.2 to 0.3 


small 5 1 ^2 



TABLE I: Allowed parameter regions in the IHDM that are consistent with the WMAP dark 
matter relic density 3a region. 



missing Et in those processes due to the undetectable S particles which are at the end of 
the decay chains. Experimentally, we can search for events with single lepton + missing 
Ej-, dilepton + missing E^, trilepton + missing E-r or in general, jets + leptons + missing 
Et- The dominant SM background comes from WW, WZ and ZZ. The collider analysis 
on this model is currently under study [37]. The scalar mass in the high mass region (V), 
however, is larger than 500 GeV. Since they only have weak interactions, the production 
cross section at the LHC is typically too small. 

For the low uih region (II), when a relatively large jA^,! ~ 0.2 could be accommodated, 
the indirect detection of the dark matter via its annihilation into neutrinos, photons, 
electrons and positrons could be very promising. A recent study on the indirect neutrino 



signals [ij] showed that in the low nis region, tens to a hundred of neutrino events per 



year from dark matter annihilation inside the Earth and hundreds of neutrino events per 
year from dark matter annihilation inside the Sun can be expected at future neutrino 
telescopes. Indirect photon signals including monochromatic photon line, fragmentation 
photori spectrum and final state radiation photon spectrum in the IHDM is under current 
study [381]. High mh regions (HI) and (IV) typically don't have promising indirect detection 
signals due to the suppression of the annihilation cross section by large m^. The indirect 
detection possibility for the high ms region (V) is also less comparing to region (II). In 
Ref. fil, it is shown that neutrino events from the Earth is too low to be observed, while 
a few events per year is expected from neutrinos from the sun. 

In summary, the IHDM is a simple extension of the SM which provides a very promising 
WIMP dark matter candidate. There are several regions of the parameter space that 
could provide the right amount of dark matter relic density in the Universe. The collider 
phenomenology of this model is very rich. For certain regions, the indirect detection via 
neutrinos or photons is possible at future neutrino and gamma ray telescopes or ground 
based experiments. 
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